Characteristics of strong ferromagnetic Josephson junctions with epitaxial barriers 
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We present the measurement of superconductor / ferromagnetic Josephson junctions, based on 
an epitaxial Nb bottom electrode and epitaxial Fe2oNiso barrier. Uniform junctions have been 
fabricated with a barrier thicknesses in the range 2 — 12 nm. The maximum critical current density 
~ 2.4 ±0.2 x 10 9 Am -2 was found for a devices with a 3 nm thick barrier at 4.2 K, corresponding 
to an average characteristic voltage IcRn ~ 16 /iV. The IcRn showed a non-monotonic behavior 
with Fe2oNi8o thickness. The variation of the resistance of a unit area ARn, of the junctions with 
barrier thickness gave a Nb/Py specific interface resistance of 6.0 ±0.5 ffim 2 and Fe2oNiso resistivity 
of 174 ± 50 nfim, consistent with other studies in polycrystalline samples. 

PACS numbers: 74.45.+C, 72.25.Mk, 75.70.Cn 
Keywords: Josephson junction, ferromagnet, epitaxial 



I. INTRODUCTION 

The early research into the proximity effect between 
superconductors (S) and ferromagnets (F) concentrated 
on measurements of the critical temperature Tq and crit- 
ical field of S/F heterostructures. A motivating factor 
for this research was the realization of the n state, in 
which the groundstate phase difference between S lay- 
ers was changes from to ir, due to the oscillation of 
the superconducting order parameter induced in the F 
layeri The transition should manifest itself as a non- 
monotonic change in the properties of the multilayers, as 
a function of F layer thickness dp. Many epitaxial and 
polycrystalline systems involving different materials were 
investigated, using various growth techniques. Although 
oscillatory Tc{(Lf) were observed in several experiments 
(see reference^ for a review) , these studies were compli- 
cated by interface effects and 'dead' magnetic layers at 
the S/F interface which made the interpretation of the 
oscillatory Tq more difficult: 3 

It was not until the measurement of current perpendic- 
ular to plane (CPP) Josephson junctions with ferromag- 
netic barriers (S/F/S) that the 7r-shift could be conclu- 
sively demonstrated. Such 7r-junctions have been char- 
acterised as a function of temperature and dp^^ using 
alloys whose composition could be tuned to achieve an 
appropriately low Curie temperature (7m), such that the 
— 7r crossover could be observed in a easily realizable 
window of experimental phase space. These 7r-junctions 
have since been incorporated into various loop geometries 
to further demonstrate the tt shift 

S/F/S junctions have also been fabricated with the rel- 
atively high Tm ferromagnets Ni^ Co^ and composite 
Co/Cu/Fe2oNi 8 o structures^ In these cases the junc- 



tions are much more sensitive to the barrier properties 
and d,F, and hence the — 7r crossover has not been 
demonstrated injunctions with strong ferromagnetic bar- 
riers. 



In all of the above CPP junctions, the S and F layers 
have been polycrystalline and in the dirty limit. The 
realization of epitaxial junctions in the clean limit may 
remove some of the difficulties of measuring low Tm alloy 
systems, (which are sensitive to stoichiometry and harder 
to characterize magnetically), as well as the sensitivity 
to dp of the high Tm barriers. Clean limit junctions 
are also expected to show a novel non-sinusoidal current- 
phase relationship^ 1 -, in contrast to recent measurements 
of S/F/S junctions with Cu^Nii-^ alloy barriers^ 



S/F/S 7r-junctions have been proposed as potential 
logic elements in quantum computing circuits 13-14 ; how- 
ever the present critical current - resistance product 
(IcRn) values are relatively small at present. The clean 
limit may again provide a route to achieving much higher 
critical current densities Jc, and hence IcRn- The com- 
bination of all of these factors motivates the investigation 
of S/F/S Josephson junctions based on epitaxial ferro- 
magnetic barriers. The epitaxial Fe2oNiso (Py) system is 
also of interest in spintronic applications, such as the fab- 
rication of ballistic spin valves and spin torque devices, 
and also provides an interesting comparison to the previ- 
ous studies of polycrystalline Nb/normal metal (N) and 
F/N interfaces (for examplei£*i£ii£) . 
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II. FILM GROWTH AND 
CHARACTERIZATION 

Our (110) Nb/(111)_ Fe 2 oNi 80 (Py) films are grown 
by sputtering on (112 0) AI2O3 as described in de- 
tail elsewhere^ To improve the epitaxy and reduce 
the strain in the films, the actual device structure is 
Nb/Cu/Py/Cu/Nb with the thicknesses of the two Cu 
layers ~ 5 nm, (with (111) orientation). The Cu layers 
are expected to be strongly proximitized by the Nb elec- 
trodes and should not significantly reduce the Jc of the 
devices. The Py thickness dp y , was in the range 2—12 
nm. The samples with dp y = 2, 4 and 6 nm were ini- 
tially grown with a bottom Nb layer of thickness 200 
nm, and a top electrode of 20 nm of Nb followed by a 
5 nm Au capping layer. This Au was removed ex-situ 
by Ar ion milling, and the top Nb electrode deposited 
by further d.c. sputtering. All of the remaining samples 
were grown in a second deposition, with the top and bot- 
tom Nb thicknesses of ~ 250 nm, deposited in-situ. To 
achieve epitaxy, the Py barrier was grown at 423 K, and 
the Nb at 1023 K. It was not possible therefore to grow 
the top Nb electrode epitaxially at the reduced tempera- 
tures required for the Cu and Py layers. Fig. ^shows the 



perconducting coherence length £5 = 18 nmji£ compared 
to a typical dirty limit value of ~ 6 nm in polycrystallinc 
sputtered filmsi 20 ! 21 
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FIG. 1: X-ray diffraction scan of (HO) Nb, (111) Cu and 
(111) Py peaks at positions 28 = 38.6, 43.2 and 43.9° respec- 
tively. The splitting of the Nb peak is caused by the presence 
of the two CuK a radiation lines. Inset: An u scan of the 
primary Nb peak. 

x-ray diffraction peaks of the Nb, Cu and Py layers in the 
rfpy = 6 nm sample, (taken with CuK Q radiation using 
a Philips X'Pert powder diffractometer) . The full width 
half maximum values obtained from u> scans were ~ 0.25° 
for the Nb layer (inset of Fig. [TJ, and ~ 0.78° for the Py 
layer (not shown). This confirms the epitaxial nature of 
the bottom electrode and barrier. A resistance vs tem- 
perature, R(T), measurement of the unpatterned dp y = 2 
nm film (before further Nb was deposited, such that the 
relatively thick epitaxial Nb layer dominates the conduc- 
tivity) was also made. The film showed a residual resis- 
tance ratio R(T = 10K)/i?(T = 300K) = 14.2. This is of 
similar order to other epitaxial Nb films, in which the su- 
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FIG. 2: Scaling of the magnetic moment per unit area vs Py 
thickness. The line is a least squares fit to the data. Inset: 
Hysteresis loop of the 2 nm Py him at 295 K. 

The magnetic properties were characterized with a vi- 
brating sample magnetometer at room temperature. Fig. 
[21 shows the saturation moment per unit area of the films 
in this study, (along with a typical hysteresis loop shown 
in the inset). Extrapolating the least square fit gives a 
nominal magnetically 'dead' layer of thickness ~ 7 A. 

The films were patterned using optical lithography, fol- 
lowed by broad beam Ar ion milling (1 mAcm~ 2 , 500 V) 
to micron scale wires and associated tracks and contact 
pads, to allow four point measurements to be performed 
on the devices. These tracks were then processed in a 
Ga + focused ion beam to achieve vertical transport with 
a device area in the range 0.05 — 1.1 /mi 2 . This fabrica- 
tion process is described in detail elsewhere 2 ^ and has 
been used previously to fabricate Josephson junctions 
with strong ferromagnetic barriers^ Transport measure- 
ments were made in a liquid He dip probe. The differen- 
tial resistance as a function of bias current of the junction 
was made with a lock- in amplifier, and the Ic found us- 
ing a resistive criterion. The Rn was measured using 
a quasi-d.c. bias current of 3 — 5 mA. This enabled the 
non-linear part of the I — V near to Ic to be neglected, 
but was not too large to drive the Nb electrodes normal. 



III. RESULTS AND DISCUSSION 

An Ic{H) modulation obtained in a junction with lat- 
eral area ~ 1060 nm x 250 nm is shown in Fig. [3] In 
this case the field is applied in the direction perpendicu- 
lar to the larger dimension of the device. It is expected 
that the coercive field of the Py should increase relative 
to that taken from the room temperature hysteresis loop 
shown in Fig. [21 due to the reduced temperature and the 
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aspect ratio of the sub-micron patterned device^ In this 
case the applied field required to modulate the Ic does 
not significantly affect the magnetization of the Py, and 
the Ic (H) is symmetric about zero field. The good fit to 
the ideal 'Fraunhofer' pattern indicates a uniform current 
flow in the junction. Relatively smaller junctions, which 
require larger fields to modulate the Ic, were found to 
show hysteresis which we associate with changes in the 
Py domain structure and magnetization. Net magnetic 
induction present in the barrier is known to shift the 
Ic(H) pattern away from H = and reduce the Iq at 
zero field from its true maximum valued To avoid falsely 
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FIG. 3: Ic(H) modulation for a device demagnetized at 4.2 K 
for a 2 nm thick Py barrier. Device dimension in the direction 
perpendicular to the applied field ~ 1.06 /im. Line is a best 
fit to a 'Fraunhofer' function. 

small values of Ic, the films were demagnetized at 4.2 
K, as well as the Ic{H) pattern being directly measured 
where possible, (for devices with IcRn > 2 /iV). In some 
cases however, the offset of the maximum Ic from H = 
could not be removed. The cause of this is ascribed to 
shape anisotropy in the sub-micron devices which pro- 
vides an additional demagnetizing field which may pre- 
vent the formation of a flux-closed domain structure in 
the barrier. 

The variation of IcRN(dp y ) is shown in Fig. 0] 
Other devices showed much larger Jc values, but showed 
strongly distorted or no Ic(H) modulation, implying 
shorting around the edges of the junctions due to rede- 
posited material during device fabrication. For dp y > 8 
nm the devices showed some reduction of the differential 
resistance around zero current bias, but did not show a 
measurable supercurrent at T = 4.2 K, (this was also 
the case in several devices with dp y — 7 and 8 nm). No 
re-entrant IcRn was observed up to dp y = 12 nm. It 
is clear that there is a strong suppression of IcRn for 
increasing dp y , but despite some scatter in the data of 
Fig. 0]for each set of devices with constant dp y , the de- 
cay is clearly not purely exponential. A component of 




4 5 

Py thickness (nm) 

FIG. 4: Characteristic voltage IcRn, as a function of Py 
thickness at T = 4.2 K. Dashed and solid lines are two fits to 
Eq. Q, as described in the text. 



this non-monotonic change may be associated with the 
slightly different preparation methods of the 2, 4 and 6 
nm thick barriers, or run to run variation in the system. 
For example for dp y = 3 nm, the fully in-situ deposition 
may imply a higher quality and larger IcRn, however 
this is inconsistent with the same comparison between 
the samples with dp y = 4, 5 and 6 nm. Therefore the 
variation in IcRn would seem to be a true effect associ- 
ated with the oscillatory induced superconducting order 
parameter in the Py layer. 

We can compare the behavior of these junctions to 
the non-epitaxial evaporated junctions of Blum et al&, 
with the structure Nb/Cu/Ni/Cu/Nb. In that case a 
measurable critical current at T = 4.2 K was observed 
up to a Ni thickness of 9 nm - similar to the present data. 
Due to the relatively small number of data points and the 
scatter, it is difficult to accurately fit this non-monotonic 
decay. We model the data using 

IcRn oc | sm{2E cx d F /hv F )\/(2E cx d F /hv F ) , (1) 

where v F is the Fermi velocity of Py, and E ex the ex- 
change energyi^ These two parameters have recently 
been measured^ in Py as E cx — 135 meV and v F — 
2.2 ±0.2 x 10 5 ms" 1 (for the majority spin). The lines in 
Fig. ^correspond to E cyi — 95 meV, with both v F = 2.2 
and 2.1 x 10 5 ms -1 , to indicate the degree of variation 
the error in v F causes. The agreement between the data 
and the model is not ideal, however for an order of mag- 
nitude estimate it is clear that the fit is acceptable. The 
period of oscillation of IcRn (d F ) is given by irhv F / E ex , 
and can therefore be estimated to be of the order of 5 nm, 
again similar in magnitude to the 5.4 nm value obtained 
by Blum et al& 

For the previously mentioned Ni junctions 8 , a max- 
imum Ic ~ 20 mA was observed at T = 4.2 K in 
10 x 10 /im devices for a 1 nm thick Ni barrier, giv- 
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ing Jc = 2 x 10 8 Am 2 . In our samples, one device with 
dp y = 3 nm had 1q = 1-03 mA with a lateral area of 
0.35 ± 0.02 fim 2 , giving the highest J c ~ 2.9 ± 0.2 x 10 9 
Am -2 at the same temperature, the average for dp y = 3 
nm was J c = 2.4 ± 0.2 x 10 9 Am" 2 . The larger J c for 
the epitaxial films can be attributed to the increased £s 
of the Nb bottom electrode. 

Finally, we have measured the total resistance (Rn) 
of a unit area ARn, of all junctions. The variation of 
ARn for all devices vs. dp y over the range of 2 nm < 
dp y < 12 nm are shown in Fig. |3J Here ARjq is the total 
specific resistances; consisting of the specific resistance 
of the S/F interfaces (AR Nh / Py ) in the S/F/S sandwich 
and the ferromagnetic layer^ such that 

ARn = 2Ai? N b/p y + ppydpy. (2) 

If we exclude one of the data points (with the highest 
ARm for dp y = 2 nm), the best fit straight line to our 
data gives an ordinate intercept of 2AR Nh / Py — 6.0 ±0.5 
iflm 2 and the resistivity pp y — 174 ± 50 nf2m (from the 
slope). The values fall within the range of 2Ai? Nb / Py = 
6 — 7.5 film 2 and pp y = 110 — 140 nflm reported for 
polycrystalline samples elsewhere iliil 
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FIG. 5: Resistance area product (ARn) at T — 4.2 K as a 
function of Py thickness. Line is a best fit excluding the point 
with the largest ARn for dp y = 2 nm. 



IV. SUMMARY 

We have measured the Josephson current through a 
epitaxial Py barrier and observed high quality junction 
characteristics. The IcRN(dp y ) showed a non-monotonic 
behavior, which could be approximately modelled with a 
simple model. We observed no sign of re-entrant behav- 
ior above dp y = 8 nm. The data extracted from the 
ARN(dpy) product was consistent with polycrystalline 
samples. 

While interesting for spin-valve junctions, Py is not 
an ideal material in which to fully explore the properties 
of ferromagnetic Josephson junctions in the clean limit 
since the spin diffusion length is relatively short A 
reduced spin diffusion length can have a strong influence 
on the possible realization of 7r-junctions formed by the 
interference of multiple Andreev reflection processes at 
the S/F interfaces. The role of different reflection ampli- 
tudes for the minority and majority spins, and the spin 
polarisation have been considered elsewhere 2&i£L Other 
effects on the Josephson current due to shape of the Fermi 
surface in the epitaxial barrier may also be a considera- 
tion, in an analogous fashion to the requirements in tun- 
nel junctions^ The use of elemental, or other epitaxial 
F layers - with longer spin diffusion lengths - may allow 
even further increases of the Jc to be achieved in these 
structures, and rule out any problems associated with a 
loss of spin memory. With a suitable choice of material 
and growth technique, the fully epitaxial S/F/S would 
also be an important system to study. 
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